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ABSTRACT 

Context. Previous studies have indicated that the 372.4 GHz ground transition of orfho-H 2 D + 
might be a powerful probe of Proto-Planetary Disks. The line could be especially suited for study 
of the disk mid-plane, where the bulk of the mass resides and where planet formation takes place. 
Aims. Provide detailed theoretical predictions for the line intensity, profile and maps expected for 
representative disk models. 

Methods. We determine the physical and chemical structure of the disks from the model devel- 
oped by Ceccarelli & Dominik (2005). The line emission is computed with the new radiative 
transfer method developed recently by Elitzur & Asensio Ramos (2006). 

Results. We present intensity maps convolved with the expected ALMA resolution, which de- 
lineate the origin of the H 2 D + 372.4 GHz line. In the disk inner regions, the line probes the 
conditions in the mid-plane out to radial distances of a few tens of AU, where Solar-like plane- 
tary systems might form. In the disk outermost regions, the line originates from slightly above 
the mid-plane. When the disk is spatially resolved, the variation of line profile across the image 
provides important information about the velocity field. Spectral profiles of the entire disk flux 
show a double peak shape at most inclination angles. 

Conclusions. Our study confirms that the 372.4 GHz H 2 D + line provides powerful diagnostics 
of the mid-plane of Proto-Planetary Disks. Current submillimeter telescopes are capable of ob- 
serving this line, though with some difficulties. The future ALMA interferometer will have the 
sensitivity to observe and even spatially resolve the H 2 D + line emission. 

Key words. ISM: abundances — ISM: molecules — stars: formation 

1. Introduction 

Proto-Planetary Disks are believed to mark the transition between the protostellar phase and planet 
formation. Similar to "standard" molecular clouds, young Proto-Planetary Disks are gas rich, with 
gas to dust mass ratio of ~100, and their dust grains have average sizes of ~0.1 yum. As the disk 
evolves with time, the dust grains coagulate into progressively larger bodies — the seeds for planet 
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Adams et al 



2004). 



formation. The gas is dispersed via viscous accretion and/or evaporation (e.g. 
Following the details of this evolution is of paramount importance for understanding the process of 
planet formation. Especially relevant in this context is the disk mid-plane because it contains the 
bulk of the mass and because its magneto-ro tational instabilities are believed to be a major source 
of the disk viscous accreti on (Gammielll996 ): thus determining the degree of ionization is also of 
great importance (see also Semenov et all2004 ). 

Observations of the 372 GHz 1 1,0 - 1 1,1 ground-state transition of ortho-H2D + provide the best, 
and possibly only, current means for studying these aspects of disks around solar type protostars. 
The reason is that in the disk mid-plane the temperature is so low (< 20 K) and the density so high 
(> 1 x 10 6 cm -3 ) that virtually all heavy-elements bearing molecules freeze-out onto the grain 
mantles, disappearing from the gas phase. Only molecules containing H and D atoms stay gaseous 
and are capable of probing the gas in the disk mid-plane. Because the ground transitions of H2 
and HD, the most abundant H and D bearing molecules, have relatively high energies (> 100 K), 
they do not get excited in that cold environment. In principle, these transitions could be seen in 
absorption against the dust continuum, but such observations are extremely difficult because, if for 
no other reason, the relevant wavelengths are blocked by the Earth atmosphere (28 and 112 pm 
for H2 and HD, respectively) Q. As a result, the most promising probes of the disk mid-plane are 
transitions of the next most abundant molecules containing H and D atoms: H3, formed by the 
interaction of cosmic rays with H2 and H, and its deuterated forms H2D + , HDJ and . While the 
symmetric molecules H^ and do not possess dipole rotational transitions, both H2D + and HDj 
have ground transitions in the sub-millimeter region that can, in principle, be observed with ground 
based telescopes. Moreover, the disappearance of heavy-element bearing molecule s from the gas 



phase enhances the H2DVH3 and HDJ/HJ ratios to values that can exceed unity dRoberts et al 



2003) — H2D + and HDt become the most abundant positive charge carriers under these circum- 



stance, and thus are not only the best probes of the disk mid-plane but also of the ionization degree. 
Detailed models of t he chemistry of Proto-Planetary Disks confi rm and quantitatively support these 



general arguments dCeccarelli & Dominikll2005 



Wil lacvll200 7). 



Both H2D + and HDj come in ortho and para forms, depending on the spin alignment of the H 
and D pairs. Of the five possible ground state transitions of t he four species, o nly three have been 



observed thus far: the ortho-H 2 D + 1 1, - In at 372.4 GHz dStark et al 



199% and the 2 L2 - 1] 



at 23 63 GHz (Cernicharo et al. 2007), and the para-HD+ 1 U) - l ,i at 691.7 GHz (IVastel et al 



2004). The two other transitions lie in a frequency domain inaccessible from the ground (at 1370. 1 



and 1476.6 GHz) and have not been detected thus far; this should change with the anticipated 
launch of Herschel, followed by SOFIA. Presently, the para-HDJ line has been detected in only 
one Pre-Stellar Core (Vastel et al. 2004). In contrast, the ortho-H2D + line at 372 GHz has been 
observed in several obj ects, mostly Pre-Stellar Cores, where the conditions are sim ilar to those 



in the disk mid-plane (Caselli et al 



2003 



Hogerheiide et al 



2006; 



Hariu et al 



2006), and in one 



Proto-Plane tary Disk (Ceccarelli et al. 2004). In Pre-Stellar Cores, th e line is so bright th at both the 



line profile dvan der Tak et al.l l2005) and the emission spatial extent (IVastel et al 



2006) have been 



studied. Unfortunately, the compact dimensions of Proto-Planetary Disks preclude such studies 
there because of the limited sensitivity and spatial resolution of current instruments. However, the 



1 Previous detections o f the H 2 S(0) and S(l) transitions towa rds Proto-Planetary Disks by 
have not been confirmed 1 Sako et al]|2005 ; Pascucci et alj|2006 ). 



Thi et al 



(2001) 
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future Atacama Large Millimeter Array (ALMA) will have the capacity to not only easily detect 
the H2D + line in disks but also to carry out mapping and line profile measurements, similar to 
what has already been done toward Pre-Stellar Cores. These studies will provide us with maps of 
the ionization structure in the disk mid-plane, and the dynamics; possibly, the presence of forming 
planets might be detected in perturbations of Keplerian motions. 

In anticipation of this progress in observational capabilities, here we report a study of the ortho- 
H2D + 1 10 - 1 1 , i line emission in Proto-Planetary Disks. We present theoretical predictions for the 
line profiles and the line emission maps in a variety of Proto-Planetary Disk models. The calcu- 
lations build on the study of the chemistry in the outer mid-plane of Proto-Planetary Disks by 



Ceccarelli & Dominik (2005), and employ the new radiative transfer formalism recently developed 
by Elitzur & Asensio Ramosl d2006l) to compute the H2D + line emission. Details of the modeling 
procedure, including non LTE radiative transfer calculations and chemical/physical models, are de- 
scribed in ^2] The results of our computations and implications for observability of the line and its 
diagnostic capability are described in |3] Our conclusions are contained in ^4]. 



2. The model 



2.1. The disk model 



Following the strategy described in lCeccarelli & Dominikl d2005l) . we use the grid of models com- 
puted in that study. We consider a Proto-Planetary Disk in Keplerian rotation around a protostar 
with mass 0.5 M Q , luminosity 0.5 L Q and T e s = 3630 K. Starting at an inner radius of 45 AU, the 



disk surface density follows a power law E 



namely, the mass per unit radius is constant 



across the disk. The gas is mixed with the dust, and we fix the mass of the disk dust content at 
2 x 10~ 4 M Q . Model number 1 corresponds to a set of "standard" parameters: disk outer radius of 
400 AU; dust-to-gas mass ratio of 1:100, i.e., the disk mass is 0.02 M ; average grain radius of 



0.1 yum; and cosmic ray ionization rate £ cl = 3 x 10 _1/ s ^. We consider four additional models, 
varying these parameters one at a time in each case. In model 2 the grain radius is increased to 1 
/mi, corresponding to a case in which a substantial fraction of the dust has coagulated so that the 
small grains are removed. Model 3 checks the effect of increasing the cosmic ray ionization rate 
by factor ten, as might be the case of either a larger f or when emulating strong X-ray irradiation 
of the disk. In model 4 the dust-to-gas ratio is increased by factor ten, so that the disk mass is only 
2xl0~ 3 M G . This corresponds to a phase in which a fair amount of gas has been dispersed. Finally, 
model 5 corresponds to a larger disk radius of 800 AU, similar to the well studied disk of DM Tau, 



where H2D + has been detected dCeccarelli et al 
models. 



2004). Table 1 summarizes the parameters of all 



Using the approach described in 



Dullemond et al 



020011) and lDuUemond & Dom inik (2004) 



for a passively irradiated hydrostatic flaring structure, the disk physical structure is derived self- 
consistently. The gas temperature is assumed to be the same as the dust. This is a justified 
assumption for the regions relevant to the present study, where the density is larger than 
about 10 5 cm 3 and the gas and dust are thermally coupled because of the gas-dust collisions 
(e.g. Goldsmith 2001). We also assume a constant dust-to-gas ratio throughout the disk, as- 



2 We assume a constant cosmic rays ionization rate across the disk, as they can penetrate up to a column 
density of ~ 100 gr cirT 2 , never exceeded in the disk regions relevant to the present study.. 



4 



Asensio Ramos et al.: H2D + line profiles 



Model 


dust: gas 


<Tcr 


^grain 


disk radius 




mass ratio 


[3 x 10- 17 s" 1 ] 


[/an] 


[AU] 


1 


0.01 


1 


0.1 


400 


2 


0.01 
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0.01 
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400 
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0.1 
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0.1 


400 
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0.01 


1 


0.1 


800 



Table 1. The five disk models considered in this study; f cr is the cosmic rays ionization rate, a, 
is the average grain radius. 



suming that no sedimentation of dust occurred yet. Fig.[TJshows the radial and vertical profiles 
of the gas density and the dust temperature for the 5 models considered here. The chemical com- 
position, in particular the H2D + radial and vertical abundance profi les, is computed by solving for 



chemical equilibrium, as described in lCeccarelli & Dominikl d2005l) . and is shown in Fig.Q]for the 



5 models. The H2D + abundance depends primarily on: i) the cosmic rays ionization rate, which 
governs the overall ionization structure; ii) the dust-to-gas ratio, which determines the disk mass; 
and iii) the grain average radius a gra i n . The grain size affects the overall ionization structure since 
recombination with negatively charged grains can become the major destruction channel for H3 , 
H2D + etc. In addition, a gra i n affects also the deuteration ratio (i.e. H2DVH3) since it controls the 
surface area on which heavy-elements bearing molecules freeze-out. The models considered here 
ex plore the influence of each re levant parameter on the H2D + li - 1q,i line. Based on the work 



bv lCeccarelli & Dom inik (2005), this ensemble of models should provide a good representation of 
the range of physical conditions relevant to H2D + line emission in Proto-Planetary Disks. 

The H2D + density, which is controlled primarily by the depletion of heavy-elements bearing 
molecules due to freeze-out, is expected to peak at "intermediate" depletion — at low depletion 
H3 becomes more abundant than H2D + , and at very large depletion D3 becomes the most abundant 
positive charge carrier, as it is the end of the chain transferring D atoms from HD. Indeed, the 
"standard" Model 1, displayed in the top row of figureQ] shows that the H2D + density is largest in 
the disk mid-plane, where the depletion is "intermediate", as expected. The "intermediate deple- 
tion" region shrinks when the grain size increases (Model 2, second row in the figure) and expands 
when the cosmic ionization rate increases (Model 3, third row); furthermore, the peak H2D + den- 
sity is higher in the latter case. Increasing the dust-to-gas ratio implies a lower gas density (as the 
disk dust mass is fixed), and a comparatively reduced region of H2D + at a given density (Model 
4). Finally, keeping the disk mass fixed while increasing its radius (Model 5) reduces the overall 
density and decreases the portion of the disk where H2D + is the most abundant; most of the outer 
disk mid-plane is now so cold that the major charge carrier there is D3 . 

2.2. The radiative transfer calculation 

At every point in the disk, the populations of the two levels of the l^o - li,i transition are de- 
termined from the balance between excitations and de-excitations due to collisional and radiative 
interactions. Given the population distribution, the intensity emerging along any line of sight can 
be obtained from ray tracing. However, the radiative interactions involve the local radiation field, 
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Fig. 1. Physical and chemical structure of the models listed in Table 1; each row corresponds to 
the model number listed in its right panel. Each panel presents a contour map whose axes denote 
radial and vertical distances normalized to the disk radius (note that Ra is doubled in Model 5). Left 
column: H2 density. Center column: Dust temperature. Right column: H2D + density. Contour units 
are cirT 3 for the H2 and H2D + densities, and K for the dust temperature. The contour values are 
denoted in the bar on top of each column. Note that Model 2 is truncated at r/R^ = 0.2, where the 
CO depletion is larger than 3; the Ceccarelli & Dominik (2005) model is no longer adequate in this 
regime. 



which includes the contribution of line emission from the entire disk (the "diffuse radiation") and 
thus cannot be determined before the level populations are known everywhere. Determining the 
level populations requires coupling between the radiative transfer and statistical rate equations. In 
the recently developed Coupled Escape Probability (CEP) method, the effects of the diffuse ra- 
diation are incorporated into the level population equation s through coefficients that rigor ously 
account for the radiative coupling between different regions ( Elitzur & Asensio Ramosll2006 ). The 
exact solution of the problem is obtained from a set of non-linear algebraic equations for the level 
populations in each zone. In a comparative study of a number of standard problems, the CEP 
method outperformed the leading accelerated A-iteration by substantial margins. 

With a grid of 80 radial points and 100 vertical points we divide the disk into 8000 rings. Each 
ring has a rectangular cross-section and uniform physical conditions. Thus far the CEP method has 
been formulated only for the slab geometry, therefore its application here involves an approxima- 
tion. Consider the vertical profiles of density, etc., in the disk at a given radius, and a stratified slab 
with the same structure perpendicular to its surface. Neglecting radiative energy flow in the radial 
direction through the disk, the vertical profile of the level population distribution at that radius can 
be obtained from the solution of the stratified slab. This is the approach we take here. The approx- 
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imation is justified because radiative coupling is suppressed in the radial direction not only by the 
large optical depths but also by the Doppler shifts due to differential Keplerian rotation. 

Following the suggestion by Stark et al. (1999) that the collisional de-excitation rate coefficient 
is ~ 1CT 10 cm 3 s _1 , the transition critical density used here is 1 x 10 6 cm" 3 . New unpublished 
computations (E. Hugo & S. Schlemmer, private communication) suggest a rate coefficient of ~ 
10~ 9 cm 3 s _1 at 10 K; changing the temperature to 20 K has only a 5% effect on the result in 
these computations. The implied critical density is a factor 10 lower than assumed here. Since 
the densities of the emitting gas are larger than either value, the transition level populations are 
practically always in thermal equilibrium so that the exact value of the collision rate is not critical. 
Finally, we assume a constant ortho-to-para H 2 D + ratio across the disk equal to 0.3, as in 
Ceccarelli & Dominik (2005). 



3. Results 

The aim of this study is to provide predictions for comparison with actual observations of the 
ortho-H2D + ground line. Two categories of measurements exist — single-dish and interferometric 
observations. We present the results of our calculations for these two cases separately in each of 
the fol lowing subsections. In compu ting the line emission, we assume the H2D + ortho-to-para ratio 
is 0.3 ( Ceccarelli & Dominikll2005 ) and consider a source at a distance of 140 pc. 



3.1. Spectral line profiles 

At present, few single-dish telescopes are capable of observing the H2D + 372 GHz line: the 15m 
JCMT, 12m APEX, 10m CSO and 3m KOSMA. At a distance of 140 pc, a disk radius of 400 
AU corresponds to an angular diameter of ~6", less than the smallest telescope beam currently 
available (JCMT with 15" at 372 GHz). The only measurable quantity is the overall flux emerging 
from the disk. Figure [2] shows the flux spectral shape for the five models listed in Table 1 for 
different viewing angles, from face-on to edge-on. The flux level is highest in Model 3, which has 
the highest H2D + density (see figure [TJ owing to its increased cosmic ray ionization rate, and in 
model 5, which has the largest surface area. It is lowest in Models 2 and 4, where the H2D + density 
is reduced. 

All models show a similar behavior for the line profiles. Face-on viewing produces single-peak 
profiles with the thermal line width. As the disk orientation changes toward edge-on, the profiles 
broaden and switch to a double-peak shape reflecting the disk rotation. The peak separation is 
determined by the combined effects of column density and velocity projection along the line of 
sight. For the 400 AU disks, the Keplerian velocity varies from 3.2 kms -1 at 45 AU to 1.1 kms" 1 
at 400 AU, with the largest line optical depth in the vertical direction (with a value of ~0.5 for 
Model 1) occurring at a radius of ~ 220 AU, where the velocity is 1.4 kms -1 . To demonstrate the 
effect of rotation speed on the profile shape, figure [2] shows also profiles in which the disk rotation 
is arbitrarily speeded up or slowed down by factor 2 while all other model parameters are held 
fixed. Although the profiles change their shape with inclination, the area under each curve remains 
almost the same; that is, the velocity-integrated flux is roughly independent of viewing angle. The 
computed values of the velo city-integrated flux agree rather well with the estimates obtained by 



Ceccarelli & Domin ik (2005) with the standard escape probability approximation. 
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Fig. 2. Spectral line profiles of the disk flux at a distance of 140 pc for each model at various 
viewing angles, as marked (face-on is 90°). For demonstrating the effects of rotation, the bottom- 
right panel reproduces the profile of Model 1 at 30° inclination together with the same calculation 
in which the Keplerian rotation is speeded up and down by factor 2, as marked. 



3.2. Line optical depth 

It is of interest to investigate the line optical depth for different models and different positions. 
Figure[3]shows the optical depth at the core of the line for models 1 and 3 for cuts passing through 
the disk center in the two limiting orientations of face-on and edge-on. Thanks to an increased 
cosmic ray ionization rate that leads to a higher H2D + density, the results of model 3 represent an 
upper limit on the optical depths. In edge-on orientation the inner 100 AU region is optically thick 
while the face-on case is neither optically thick nor thin, except for a ring at ~250 AU. This ring can 
be recognized in the images shown in Fig. |4j even though the inclination angle is different. Model 
1 produces very similar results but with reduced optical depths. In this case, only the edge-on 
orientation yields optical depths above unity. 



3.3. Brightness maps 

Interferometic observations with SMA (and eSMA), IRAM Plateau de Bure and ALMA will be 
able to spatially resolve the H2D + 372 GHz line emission from disks like the ones considered 
here. We produce model maps by convolving the intensity with a Gaussian beam with FWHM 
of 0.3", corresponding to a linear size of 40 AU at a distance of 140 pc; such resolution can be 
reasonably expected at ALMA for the predicted intensities. Fig. [4] shows some results. The edge- 
on map of Model 1 traces closely the horizontal and vertical distributions of the H2D + density. The 
most intense emission is produced in the inner regions close to the disk mid-plane within a height 
of ~ 1/10 of the radius, confirming the suggestion that the ortho-H2D + 372 GHz line probes the 
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Fig. 3. Optical depths at the core of the line for face-on and edge-on orientations. Shown are cuts 
along axes passing through the center of the disk, so that the optical depth is maximized. The results 
for model 3 represent an upper limit for the optical depth among all models. Model 1 presents 
optical depths that are intermediate between model 3 and all other models. Note that the zero 
optical depth in the central region of the face-on disk is a numerical artifact, as our model 
does not consider that region. 



disk mid-plane dCeccarelli et al.ll2004l) . The high brightness is maintained out to radii approaching 
100 AU, covering the region where a solar-like planetary system might be forming. Therefore this 
line can serve as a probe of the ionization degree and kinematics in a region critical for planet 
formation. Further out, because of the extreme molecular depletion in the equatorial plane, which 
shifts all deuterium into (see Fig.[T]and the discussion in 32. Il l, the H2D + line probes regions 
slightly above it. 

The five model maps at 30° inclination show additional structure, reflecting variation in column 
density and path length due to the geometry and chemistry of the flaring disk. These variations 
control also the channel maps, shown in Fig.|5]for the standard disk (Model 1) at 30° inclination. 
Brighter and fainter rings are evident at different positions for different velocities, depending on 
the intercepted portion of disk. Overall, the maps show the standard Keplerian pattern with blue- 
shifted and red-shifted lobes spatially separated and symmetric with respect to the rotation axis. 
Fig.|6]shows a mosaic of the line profiles for this case. 

All maps employ intensity units. Those can be converted to "equivalent Rayleigh- 
Jeans temperature" through r RJ = 2.35/ 13 K, where / n is the intensity in units of 
10 13 erg cm 2 s 1 Hz 1 sr It should be noted that, since the temperature equivalent of 
the 372 GHz transition frequency (806 //m wavelength) is 18 K, the actual brightness temper- 
ature is usually quite different from T RJ . 
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Fig. 4. Maps averaged over a representative ALMA beam (Gaussian with FWHM = 0.3") of the 
velocity-integrated intensity. The bar to the right of each panel shows the contour scales. The x and 
y axes are displacements, in AU, from the center of the disk. Top left panel is for an edge-on (90°) 
view of Model 1. The other panels show all five disk models (see Table 1) at a viewing angle of 
30°. Note the larger spatial extent of Model 5. 



4. Conclusions 



Our results confirm the utility of the 372 GHz H2D + line in probing the inner regions of Proto- 
Planetary Disks. While the predicted line fluxes are at the limit of detectable power of current 
sub-millimeter telescopes, the future ALMA interferometer will be able to not only detect the line 
but also study the distribution and extent of the H2D + emission. The spatially and velocity inte- 
grated intensities computed i n this study (§|3} compare rather well with the previous estimates by 
Ceccarelli & Dominik (2005), which were derived in the escape probability approximation. Indeed, 
Elitzur & Asensio Ramosl J2006I) note that the escape probability approximation is adequate for cal- 
culating the line integrated emission of a two-level system from a homogeneous slab. Our study 
suggests that this approximation can be used to provide rough estimates of the integrated line in- 
tensity also in the case of disks with variable physical properties. However, line profile calculations 
require an exact formalism, like the CEP method employed here. 

The major unc ertainty in our modeling is t he disk physical structure and H2D + abundance. 
As pointed out bv lCeccarelli & Dominikl (120051) . one major uncertainty is linked to the role of N2 
and, specifically, its freezing-out onto the grain mantles. The models studied here adopt a sticking 
coeffi cient of 1 and a b inding energy of 575 K, as suggested by early models and observations 
(e.g. iBergin et al.ll2002l) . However, laboratory studies show that the binding energy of N2 is sim- 
ilar to that of CO, namely 885 K. If this is the case, N2 will freeze-out more easily and a larger 
region of the disk than computed here would be dominated by DJ, at the expense of H2D + . The 
resulting H2D + column density would then be a factor 30 lower for th e standard case (Model 1), 
making the line much more difficult to detect (see also IWillacvll2007h . Complementary observa- 
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Fig. 5. Channel maps of the H2D + line intensity (convolved with ALMA beam; see Fig. |4j for the 
standard disk (Model 1), viewed at an inclination of 30°. The x and y axes are displacements, in 
AU, from the disk center. The central velocity of each map is marked on top. The bar to the right 
of each panel shows the contour brightness scale. Note the scale changes among the panels. 



tions of N-bearing molecules, like N2H + , will be necessary for drawing a coherent picture. A sec- 
ond uncertainty involves the H2D + ortho-to-para ratio, which is observationally entirely unknown. 



ortho -to-para ratio of H2 and the grain size (Flower et al 



(Flower et al. 


2004; 


Vastel et al. 


2006; 


Flower et al. 



2006). Clarifying this issue requires observations of the ground transition of the para-H2D + at 



1370.1 GHz. Unfortunately, no current facility is capable of such observations, and neither will 
Herschel. A third uncertainty arises from the collision rate coefficient of the 1 - lio transition, 
but this should have only a minor effect given the relatively high densities in the disk mid-plane. 

With all of these caveats and in spite of the uncertainties, observations of H2D + at 372 GHz 
remain a powerful diagnostic tool for probing the disk mid-plane. Imaging of the 372 GHz line 
emission will bring information, otherwise unavailable, on the disk kinematics and ionization struc- 
ture, serving as a unique tool in the study of planet formation. For example, ALMA will be able 
to reach a rms of ~2 mJy on a 0.3 km/s velocity bin and 0.3" beam in 10 hr of integration timj^, 
corresponding to a line brightness of 1 X 10 14 erg cm~ 2 s _1 Hz 1 sr~' . Comparison with the predic- 



3 based on the ALMA time estimator; http://www.eso.org/projects/alma/science/bin/sensitivity.html 
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Fig. 6. Variation of line profile across the disk image for Model 1 at various inclinations, as marked. 
The displayed mosaic is centered on the disk image and was produced as follows: The maps shown 
in figure[4]are divided into grids of 5x5 squares, each with 1.2" on the side. At a given velocity, 
the brightness is averaged over each square with a Gaussian beam of 0.3", the ALMA beam size. 



tions reported in Fig. [5] shows that this will allow imaging of the H2D + line emission at the scale 
expected for planet formation, providing a crucial information for theories of planetary formation. 
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ABSTRACT 



Context. Previous studies have indicated that the 372.4 GHz ground transition of ortho-H2D + might be a powerful probe of Proto- 
Planetary Disks. The line could be especially suited for study of the disk mid-plane, where the bulk of the mass resides and where 
planet formation takes place. 

Aims. Provide detailed theoretical predictions for the line intensity, profile and maps expected for representative disk models. 
Methods. We determine the physical and chemical structure of the disks from the model developed by Ceccarelli & Dominik (2005). 
The line emission is computed with the new radiative transfer method developed recently by Elitzur & Asensio Ramos (2006). 
Results. We present intensity maps convolved with the expected ALMA resolution, which delineate the origin of the H 2 D + 372.4 
GHz line. In the disk inner regions, the line probes the conditions in the mid-plane out to radial distances of a few tens of AU, where 
Solar-like planetary systems might form. In the disk outermost regions, the line originates from slightly above the mid-plane. When 
the disk is spatially resolved, the variation of line profile across the image provides important information about the velocity field. 
Spectral profiles of the entire disk flux show a double peak shape at most inclination angles. 

Conclusions. Our study confirms that the 372.4 GHz H2D + line provides powerful diagnostics of the mid-plane of Proto-Planetary 
Disks. Current submillimeter telescopes are capable of observing this line, though with some difficulties. The future ALMA interfer- 
ometer will have the sensitivity to observe and even spatially resolve the H 2 D + line emission. 

Key words. ISM: abundances — ISM: molecules — stars: formation 



1. Introduction 

Proto-Planetary Disks are believed to mark the transition be- 
tween the protostellar phase and planet formation. Similar to 
"standard" molecular clouds, young Proto-Planetary Disks are 
gas rich, with a gas to dust mass ratio of -100, and their dust 
grains have average sizes of ~0. 1 fim. As the disk evolves with 
time, the dust grains coagulate into progressively larger bod- 
ies — the seeds for planet formation. The gas is dispersed vi a 
viscous accretion and/or evaporation (e.g. lAdams et al.l [2004). 
Following the details of this evolution is of paramount im- 
portance for understanding the process of planet formation. 
Especially relevant in this context is the disk mid-plane be- 
cause it contains the bulk of the mass and because its magneto- 
rotational instabilities are believed to be a major source of 
the disk viscous accretion (Gammie 1996); thus determining 
the degree of ioniza tion is also of great importance (see also 
ISemenov et ai]|2004l) . 

Observations of the 372 GHz li o - li,i ground-state transi- 
tion of ortho-H2D + provide the best, and possibly only, current 
means for studying these aspects of disks around solar type pro- 
tostars. The reason is that in the disk mid-plane the temperature 
is so low (< 20 K) and the density so high (> 1 x 10 6 cm -3 ) 
that virtually all heavy-elements bearing molecules freeze-out 
onto the grain mantles, disappearing from the gas phase. Only 
molecules containing H and D atoms stay gaseous and are capa- 
ble of probing the gas in the disk mid-plane. Because the ground 
transitions of H2 and HD, the most abundant H and D bearing 



molecules, have relatively high energies (> 100 K), they do not 
get excited in that cold environment. In principle, these transi- 
tions could be seen in absorption against the dust continuum, 
but such observations are extremely difficult because, if for no 
other reason, the relevant wavelengths are blocked by the Earth 
atmosphere (28 and 112 /mi for H2 and HD, respectively) Q. As 
a result, the most promising probes of the disk mid-plane are 
transitions of the next most abundant molecules containing H 
and D atoms: H3 , formed by the interaction of cosmic rays with 
H 2 and H, and its deuterated forms H 2 D + , HDJ and D3 . While 
the symmetric molecules H^ and D^ do not possess dipole rota- 
tional transitions, both H 2 D + and HDj have ground transitions 
in the sub-millimeter region that can, in principle, be observed 
with ground-based telescopes. Moreover, the disappearance of 
heavy-element bearing molecules from the gas phase enhances 
the H2D + /H:t and HD ^/H, ratios to values that can exceed unity 
dRobertsetal.1 120031) . H 2 D + and HD+ become the most abun- 
dant positive charge carriers under these circumstance, and thus 
are not only the best probes of the disk mid-plane but also of 
the ionization degree. Detailed models of the chemistry of Proto- 
Planetary Disks confirm and quantitatively support the se general 
arguments dCeccarelli & Dom inik 2005; Willacv 2003). 

Both H 2 D + and HDj come in ortho and para forms, depend- 
ing on the spin alignment of the H and D pairs. Of the five pos- 
sible ground state transitions of the four species, only three have 
been observed thus far: the ortho-H 2 D + 1 10 - 1 1,1 at 372.4 GHz 
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1 Previous detections o f the H 2 S (0) an d S(l) transitions towards 
Proto-Planetary Disks by | Thi et alj 2001) have not been confirmed 
iSako et alJIfoollPascucci et al.ll2 006). 



(IStarket al.ll 19991) and the 2 U - 1 u at 2363 GHz (C ernicharo et 
al. 20 07), and the para-HD+ ii,o-l'o,i at 691.7 GHz (IVastel et al.l 



2004). The two other transitions lie in a frequency domain inac- 
cessible from the ground (at 1370.1 and 1476.6 GHz) and have 
not been detected thus far; this should change with the antic- 
ipated launch of Herschel, followed by SOFIA. Presently, the 
para-HDj line has been detected in only one Pre-Stellar Core 
(Vastel et al. 2004). In contrast, the ortho-H 2 D + line at 372 
GHz has been observed in several objects, mostly Pre-Stellar 
Cores, where the conditions are similar to those in the disk mid- 
plane dCaselli et all 120031; iHogerheiide et al 1 120061; lHarju et all 
120061) . and in one Proto-Planetary Disk (ICeccarelli et al.ll2004l) . 
In Pre- Stellar Cores, the line is so bright that both the line 
profile ( van der Tak et al.|[2005l) and the emission spatial extent 
( Vast el et al. 2006) have been studied. Unfortunately, the com- 
pact dimensions of Proto-Planetary Disks preclude such studies 
because of the limited sensitivity and spatial resolution of cur- 
rent instruments. However, the future Atacama Large Millimeter 
Array (ALMA) will have the capacity to not only easily detect 
the H 2 D + line in disks but also to carry out mapping and line 
profile measurements, similar to what has already been done to- 
ward Pre-Stellar Cores. These studies will provide us with maps 
of the ionization structure in the disk mid-plane, and the dynam- 
ics; possibly, the presence of forming planets might be detected 
in perturbations of Keplerian motions. 

In anticipation of this progress in observational capabil- 
ities, we report here a study of the ortho-H2D + lj o - lj i 
line emission in Proto-Planetary Disks. We present theoreti- 
cal predictions for the line profiles and the line emission maps 
in a variety of Proto-Planetary Disk models. The calculations 
build on the study of th e chemistry in the outer mid -plane of 
Proto-Planetary Disks bv lCeccarelli & Dominikl d2005l) . and em- 
pl oy the new radiative trans fer formalism recently developed 
by lElitzur & Asensio Ramosl d2006l) to compute the H 2 D + line 
emission. Details of the modeling procedure, including non LTE 
radiative transfer calculations and chemical/physical models, are 
described in Sec. [2] The results of our computations and impli- 
cations for observability of the line and its diagnostic capability 
are described in Sec. [3] Our conclusions are contained in Sec. [4] 



2. The model 

2. 1 . The disk model 



Follow ing the strategy described in ICeccarelli & Dom inik 
J20051) . we use the grid of models computed in that study. We 
consider a Proto-Planetary Disk in Keplerian rotation around a 
protostar with mass 0.5 M , luminosity 0.5 L and T e ff = 3630 
K. Starting at an inner radius of 45 AU, the disk surface density 
follows a power law E oc r _1 , namely, the mass per unit radius is 
constant across the disk. The gas is mixed with the dust, and we 
fix the mass of the disk dust content at 2 x 10~ 4 M Q . Model num- 
ber 1 corresponds to a set of "standard" parameters: disk outer 
radius of 400 AU; dust-to-gas mass ratio of 1:100, i.e., the disk 
mass is 0.02 M ; average grain radius of 0.1 /mi; and cosmic 
ray ionization rate f cl = 3 x 10~ 17 s _1 0. We consider four ad- 
ditional models, varying these parameters one at a time in each 
case. In model 2 the grain radius is increased to 1 /mi, corre- 
sponding to a case in which a substantial fraction of the dust has 
coagulated so that the small grains are removed. Model 3 checks 



the effect of increasing the cosmic ray ionization rate by factor 
ten, as might be the case of either a larger £ or when emulating 
strong X-ray irradiation of the disk. In model 4 the dust-to-gas 
ratio is increased by factor ten, so that the disk mass is only 
2xl0~ 3 M . This corresponds to a phase in which a lot of gas 
has been dispersed. Finally, model 5 corresponds to a larger disk 
radius of 800 AU, similar to th e well-studied disk of DM Tau, 
where H 2 D + has been detected dCeccarelli et alj|2004l) . Table 1 
summarizes the parameters of all models. 



Model 


dust:gas 






disk radius 




mass ratio 


[3 x 10~ 17 s-»] 


[jum] 


[AU] 


1 


0.01 


1 


0.1 


400 


2 


0.01 


1 


1.0 


400 


3 


0.01 


10 


0.1 


400 


4 


0.1 


1 


0.1 


400 


5 


0.01 


1 


0.1 


800 



Table 1. The five disk models considered in this study; £ cl is the 
cosmic rays ionization rate, a gra i n is the average grain radius. 



Using the approach describ ed in lDullemond et all ([200 1 ) and 
Dullemond & Dominikl d2004) for a passively irradiated hydro- 
static flaring structure, the disk physical structure is derived self- 
consistently. The gas temperature is assumed to be the same as 
the dust. This is a justified assumption for the regions relevant 
to the present study, where the density is larger than about 10 5 
crrT 3 and the gas and dust are thermally coupled because of 
the gas-dust collisions (e.g. Goldsmith 2001). We also assume 
a constant dust-to-gas ratio throughout the disk, assuming that 
no sedimentation of dust occurred. Figure Q] shows the radial 
and vertical profiles of the gas density and the dust tempera- 
ture for the 5 models considered here. The chemical composi- 
tion, in particular the H 2 D + radial and vertical abundance pro- 
files, is c omputed by solving for chem ical equilibrium, as de- 
scribed in ICeccarelli & Dominikl d2005l) . and is shown in Fig. Q] 
for the 5 models. The H 2 D + abundance depends primarily on: 
i) the cosmic rays ionization rate, which governs the overall 
ionization structure; ii) the dust-to-gas ratio, which determines 
the disk mass; and iii) the grain average radius a gra i n . The grain 
size affects the overall ionization structure since recombination 
with negatively charged grains can become the major destruc- 



tion channel for Hi 



H,D + etc. 
H 2 D + 



2 We assume a constant cosmic rays ionization rate across the disk, 
as they can penetrate up to a column density of ~ 100 gr cirT 2 , never 
exceeded in the disk regions relevant to the present study. 



i3 , ii 2 i 

the deuteration ratio (i.e. H 2 D + /Hp since it controls the sur- 
face area on which heavy-elements bearing molecules freeze- 
out. The models considered here explore the influence of each 
relevant parameter on the H 2 D + Up - lo. i line. Based on the 
work by ICeccarelli & Dominikl (|2005), this ensemble of mod- 
els should provide a good representation of the range of physi- 
cal conditions relevant to H 2 D + line emission in Proto-Planetary 
Disks. 

The H 2 D + density, which is controlled primarily by the de- 
pletion of heavy-elements bearing molecules due to freeze-out, 
is expected to peak at "intermediate" depletion — at low deple- 
tion HJ becomes more abundant than H 2 D + , and at very large 
depletion D3 becomes the most abundant positive charge car- 
rier, as it is the end of the chain transferring D atoms from HD. 
Indeed, the "standard" Model 1, displayed in the top row of Fig. 
[T| shows that the H 2 D + density is largest in the disk mid-plane, 
where the depletion is "intermediate", as expected. The "inter- 
mediate depletion" region shrinks when the grain size increases 
(Model 2, second row in the figure) and expands when the cos- 
mic ionization rate increases (Model 3, third row); furthermore, 



r/K, 

Fig. 1. Physical and chemical structure of the models listed in Table 1; each row corresponds to the model number listed in its right 
panel. Each panel presents a contour map whose axes denote radial and vertical distances normalized to the disk radius (note that 
Ri is doubled in Model 5). Left column: H2 density. Center column: Dust temperature. Right column: H2D + density. Contour units 
are cirT 3 for the H2 and H2D + densities, and K for the dust temperature. The contour values are denoted in the bar on top of each 
column. Note that Model 2 is truncated at r/R^ = 0.2, where the CO depletion is larger than 3; the Ceccarelli & Dominik (2005) 
model is no longer adequate in this regime. 



the peak H2D + density is higher in the latter case. Increasing the 
dust-to-gas ratio implies a lower gas density (as the disk dust 
mass is fixed), and a comparatively reduced region of H2D + at 
a given density (Model 4). Finally, keeping the disk mass fixed 
while increasing its radius (Model 5) reduces the overall density 
and decreases the portion of the disk where H2D + is the most 
abundant; most of the outer disk mid-plane is now so cold that 
the major charge carrier there is D3 . 

2.2. The radiative transfer calculation 

At every point in the disk, the populations of the two levels of 
the lio - 1 14 transition are determined from the balance be- 
tween excitations and de-excitations due to collisional and ra- 
diative interactions. Given the population distribution, the in- 
tensity emerging along any line of sight can be obtained from 
ray tracing. However, the radiative interactions involve the lo- 
cal radiation field, which includes the contribution of line emis- 
sion from the entire disk (the "diffuse radiation") and thus can- 
not be determined before the level populations are known ev- 
erywhere. Determining the level populations requires coupling 
between the radiative transfer and statistical rate equations. 
In the recently developed Coupled Escape Probability (CEP) 
method, the effects of the diffuse radiation are incorporated into 
the level population equations through coefficients that rigor- 
ously account for the radiative couplin g between different re- 
gions (lElitzur & Asensio Ramosl 12006). The exact solution of 
the problem is obtained from a set of non-linear algebraic equa- 
tions for the level populations in each zone. In a comparative 
study of a number of standard problems, the CEP method outper- 
formed the leading accelerated A-iteration by substantial mar- 
gins. 

With a grid of 80 radial points and 100 vertical points we 
divide the disk into 8000 rings. Each ring has a rectangular 
cross-section and uniform physical conditions. Thus far the CEP 
method has been formulated only for the slab geometry, there- 
fore its application here involves an approximation. Consider the 
vertical profiles of density, etc., in the disk at a given radius, 
and a stratified slab with the same structure perpendicular to its 



surface. Neglecting radiative energy flow in the radial direction 
through the disk, the vertical profile of the level population dis- 
tribution at that radius can be obtained from the solution of the 
stratified slab. This is the approach we take here. The approx- 
imation is justified because radiative coupling is suppressed in 
the radial direction not only by the large optical depths but also 
by the Doppler shifts due to differential Keplerian rotation. 

Following the suggestion by Stark et al. (1999) that the colli- 
sional de-excitation rate coefficient is ~ 10~ 1() cm 3 s~', the transi- 
tion critical density used here is 1 x 10 6 cirT 3 . New unpublished 
computations (E. Hugo & S. Schlemmer, private communica- 
tion) suggest a rate coefficient of ~ 10~ 9 cm 3 s _1 at 10 K; chang- 
ing the temperature to 20 K has only a 5% effect on the result 
in these computations. The implied critical density is a factor 
10 lower than assumed here. Since the densities of the emitting 
gas are larger than either value, the transition level populations 
are practically always in thermal equilibrium so that the exact 
value of the collision rate is not critical. Finally, we assume a 
constant ortho-to-para H2D + ratio across the disk equal to 0.3, 
as in Ceccarelli & Dominik (2005). 

3. Results 

The aim of this study is to provide predictions for comparison 
with actual observations of the ortho-H2D + ground line. Two cat- 
egories of measurements exist — single-dish and interferometric 
observations. We present the results of our calculations for these 
two cases separately in each of the following subsections. In 
computing t he line emission, we assume the H2D + ortho-to-para 
ratio is 0.3 (ICeccarelli & Dominikll2005l) and consider a source 
at a distance of 140 pc. 

3.1. Spectral line profiles 

At present, few single-dish telescopes are capable of observing 
the H 2 D + 372 GHz line: the 15m JCMT, 12m APEX, 10m CSO 
and 3m KOSMA. At a distance of 140 pc, a disk radius of 400 
AU corresponds to an angular diameter of ~6", less than the 
smallest telescope beam currently available (JCMT with 15" 



Fig. 2. Spectral line profiles of the disk flux at a distance of 140 pc for each model at various viewing angles, as marked (face-on is 
90°). For demonstrating the effects of rotation, the bottom-right panel reproduces the profile of Model 1 at 30° inclination together 
with the same calculation in which the Keplerian rotation is speeded up and down by factor 2, as marked. 



at 372 GHz). The only measurable quantity is the overall flux 
emerging from the disk. Figure [2] shows the flux spectral shape 
for the five models listed in Table 1 for different viewing an- 
gles, from face-on to edge-on. The flux level is highest in Model 
3, which has the highest H2D + density (Fig.[TJ owing to its in- 
creased cosmic ray ionization rate, and in model 5, which has 
the largest surface area. It is lowest in models 2 and 4, where the 
H2D + density is reduced. 

All models show a similar behavior for the line profiles. 
Face-on viewing produces single-peak profiles with the thermal 
line width. As the disk orientation changes toward edge-on, the 
profiles broaden and switch to a double-peak shape reflecting 
the disk rotation. The peak separation is determined by the com- 
bined effects of column density and velocity projection along the 
line of sight. For the 400 AU disks, the Keplerian velocity varies 
from 3.2 kms -1 at 45 AU to 1.1 kms -1 at 400 AU, with the 
largest line optical depth in the vertical direction (with a value 
of ~0.5 for Model 1) occurring at a radius of ~ 220 AU, where 
the velocity is 1.4 kms" 1 . To demonstrate the effect of rotation 
speed on the profile shape, Fig.|2]shows also profiles in which the 
disk rotation is arbitrarily speeded up or slowed down by factor 
2 while all other model parameters are held fixed. Although the 
profiles change their shape with inclination, the area under each 
curve remains almost the same; that is, the velocity-integrated 
flux is roughly independent of viewing angle. The computed val- 
ues of the velocity-i ntegrated flux agree ra t her w ell with the es- 
timates obtained by ICeccarelli & Dominikl d2005l) with the stan- 
dard escape probability approximation. 

3.2. Line optical depth 

It is of interest to investigate the line optical depth for different 
models and different positions. Figure [3] shows the optical depth 
at the core of the line for models 1 and 3 for cuts passing through 
the disk center in the two limiting orientations of face-on and 
edge-on. Thanks to an increased cosmic ray ionization rate that 
leads to a higher H2D + density, the results of model 3 represent 
an upper limit on the optical depths. In edge-on orientation the 
inner 100 AU region is optically thick while the face-on case is 
neither optically thick nor thin, except for a ring at ~250 AU. 
This ring can be recognized in the images shown in Fig. [4] even 
though the inclination angle is different. Model 1 produces very 



similar results but with reduced optical depths. In this case, only 
the edge-on orientation yields optical depths above unity. 

3.3. Brightness maps 

Interferometic observations with SMA (and eSMA), IRAM 
Plateau de Bure and ALMA will be able to spatially resolve the 
H2D + 372 GHz line emission from disks like those considered 
here. We produce model maps by convolving the intensity with 
a Gaussian beam with FWHM of 0.3", corresponding to a lin- 
ear size of 40 AU at a distance of 140 pc; such resolution can 
be reasonably expected at ALMA for the predicted intensities. 
Fig. 2] shows some results. The edge-on map of Model 1 traces 
closely the horizontal and vertical distributions of the H2D + den- 
sity. The most intense emission is produced in the inner regions 
close to the disk mid-plane within a height of ~ 1/10 of the ra- 
dius, confirming the suggestio n that the ortho-fyD^ 372 GHz 
line probes the disk mid-plane dCeccarelli et alj |2004). The high 
brightness is maintained out to radii approaching 100 AU, cov- 
ering the region where a solar-like planetary system might be 
forming. Therefore this line can serve as a probe of the ioniza- 
tion degree and kinematics in a region critical for planet forma- 
tion. Further out, because of the extreme molecular depletion in 
the equatorial plane, which shifts all deuterium into (see Fig. 
Q]and the discussion in Sec. 12.11 ). the H2D + line probes regions 
slightly above it. 

The five model maps at 30° inclination show additional 
structure, reflecting variation in column density and path length 
due to the geometry and chemistry of the flaring disk. These 
variations also control the channel maps, shown in Fig.|5]for the 
standard disk (Model 1) at 30° inclination. Brighter and fainter 
rings are evident at different positions for different velocities, 
depending on the intercepted portion of disk. Overall, the maps 
show the standard Keplerian pattern with blue-shifted and red- 
shifted lobes spatially separated and symmetric with respect to 
the rotation axis. Fig. [6] shows a mosaic of the line profiles for 
this case. 

All maps employ intensity units. These can be con- 
verted to "equivalent Rayleigh-Jeans temperature" through 
Trj = 2.35 7_i3 K, where 7_b is the intensity in units of 
10~ 13 erg crrT 2 s _1 Hz 1 sr _1 . It should be noted that since the 
temperature equivalent of the 372 GHz transition frequency (806 




Fig. 3. Optical depths at the core of the line for face-on and edge-on orientations. Shown are cuts along axes passing through the 
center of the disk, so that the optical depth is maximized. The results for model 3 represent an upper limit for the optical depth 
among all models. Model 1 presents optical depths that are intermediate between model 3 and all other models. Note that the zero 
optical depth in the central region of the face-on disk is a numerical artifact, as our model does not consider that region. 
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Fig. 4. Maps averaged over a representative ALMA beam (Gaussian with FWHM = 0.3") of the velocity-integrated intensity. The 
bar to the right of each panel shows the contour scales. The x and y axes are displacements, in AU, from the center of the disk. Top 
left panel is for an edge-on (90°) view of Model 1. The other panels show all five disk models (see Table 1) at a viewing angle of 
30°. Note the larger spatial extent of Model 5. 



pm wavelength) is 18 K, the actual brightness temperature is 
usually quite different from Trj. 



4. Conclusions 

Our results confirm the utility of the 372 GHz H2D + line in prob- 
ing the inner regions of Proto-Planetary Disks. While the pre- 
dicted line fluxes are at the limit of detectable power of current 
sub-millimeter telescopes, the future ALMA interferometer will 
be able to not only detect the line but also study the distribu- 
tion and extent of the H2D + emission. The spatially and veloc- 
ity integrated intensities computed in th is study (Sec. [3]l com- 
pare w ell with the previous estimates by ICeccarelli & Dominik 
(2005), which were derived in the escape proba bility approxi- 
mation. Indeed, lElitzur & Asensio Ramos (2006) note that the 
escape probability approximation is adequate for calculating the 
line integrated emission of a two-level system from a homoge- 
neous slab. Our study suggests that this approximation can be 
used to provide rough estimates of the integrated line inten- 
sity also in the case of disks with variable physical properties. 



However, line profile calculations require an exact formalism, 
like the CEP method employed here. 

The major uncertainty in our modeling is the disk phys- 
ical structure and H^D" 1 " abundance. As pointed out by 
ICeccarelli & Dominikl (120051) . one major uncertainty is linked 
to the role of N2 and, specifically, its freezing-out onto the grain 
mantles. The models studied here adopt a sticking coefficient of 
1 and a binding energy of 575 K, as s uggested by early mod- 
els and observations (e.g. Bergi rTetail 12002b . However, labo- 
ratory studies show that the binding energy of N2 is similar to 
that of CO, namely 885 K. If this is the case, N2 will freeze-out 
more easily and a larger region of the disk than computed here 
would be dominated by D3 , at the expense of H2D + . The result- 
ing H2D + column density would then be a factor 30 lower for 
the standard case (model 1), making the line much more dif- 
ficult to detect (see also IWillacvll2007l) . Complementary ob- 
servations of N-bearing molecules, like N2H + , will be neces- 
sary to draw a coherent picture. A second uncertainty involves 
the H2D + ortho-to-para ratio, which is observationally entirely 
unknown. Theoretical estimates vary depending on several pa- 
rameters, most of them unknown too, such as the ortho-to-para 




Fig. 5. Channel maps of the H2D + line intensity (convolved with ALMA beam; see Fig. [4j> for the standard disk (Model 1), viewed 
at an inclination of 30°. The x and y axes are displacements, in AU, from the disk center. The central velocity of each map is marked 
on top. The bar to the right of each panel shows the contour brightness scale. Note the scale changes among the panels. 




Fig. 6. Variation of line profile across the disk image for Model 1 at various inclinations, as marked. The displayed mosaic is centered 
on the disk image and was produced as follows: The maps shown in Fig.|4]are divided into grids of 5x5 squares, each with 1.2" on 
the side. At a given velocity, the brightness is averaged over each square with a Gaussian beam of 0.3", the ALMA beam size. 



ratio of H2 and the grai n size dFlower et al.l l2004t : IVastel et alj 
2006: Flower et al.l2006l) . Clarifying this issue requires observa- 
tions of the ground transition of the para-H2D + at 1370.1 GHz. 
Unfortunately, no current facility is capable of such observa- 
tions, and neither will Herschel. A third uncertainty arises from 
the collision rate coefficient of the 1 1 ;1> - 1 j transition, but this 
should have only a minor effect given the relatively high densi- 
ties in the disk mid-plane. 



With all of these caveats and in spite of the uncertainties, 
observations of H2D + at 372 GHz remain a powerful diagnos- 
tic tool for probing the disk mid-plane. Imaging of the 372 GHz 
line emission will bring information, otherwise unavailable, on 
the disk kinematics and ionization structure, serving as a unique 
tool in the study of planet formation. For example, ALMA will 
be able to reach a rms of ~2 mJy on a 0.3 km/s velocity bin 



and 0.3" beam in 10 hr of integration tim^| corresponding to a 
line brightness of 1 x 10~ 14 erg cirT 2 s _1 Hz -1 sr _1 . Comparison 
with the predictions reported in Fig. [5] shows that this will al- 
low imaging of the H2D + line emission at the scale expected for 
planet formation, providing crucial information for theories of 
planetary formation. 
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